Behavior of cracked materials 
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Abstract. Due to their microstructure, quasi brittle materials present rough cracks. Under sliding of 
the crack lips, this roughness involves in one hand induced opening and in the other hand some ap- 
parent plasticity which is due to the interlocking of the crack lips combined with Coulomb's friction. 
The proposed model is written under the irreversible thermodynamics framework. Micromechanics 
uses the Del Piero and Owen's structured deformation theory. Opening of the crack depends upon 
the crack shape and the relative sliding of the crack lips. The thermodynamic force associated to the 
sliding has the mechanical meaning of the force acting in order to make the crack slide. Yield surface 
is defined as a limitation of this force with respect to the Coulomb's friction and the Barenblatt cohe- 
sion. The crack orientation is defined as the one for which the criterion is reached for the lowest stress 
level. A decreasing cohesion, respect to sliding is supposed. Tension and compression reference cases 
are envisaged. 

Introduction 

This work is in continuity with J3]|, in which the problem of sliding rough fracture is treated is 
the case of compressive and shearing testing, relevant to the behavior of civil engineering structure 
response under earthquake. The goal was to extend this approach to cases of tension and compression, 
predict the crack's orientation and the global response of the specimen. 

Micromechanics 

In case of periodic fracture, a micro-mechanical structured deformation analysis 01 [31 lead to a 




Figure 1 : sketch of the micromechanics 



partition of the strain e between it's elastic (reversible) e e and due to sliding (dissipative) e s parts (Eq. 
Q]). The sliding of the crack is defined by a variable s and the opening respect to the function 9(s). 



Considering the orthonormal basis [e 1; e 2 ], where e 2 defines the normal to the (average) crack plane 
and ei the sliding direction. 

(1) 



e e + e s 



(ei ® e x ) + s9(s)~ (ei 



e 2 + e 2 <g> ei . 



(2) 



The micro-mechanical point of view is sketched in the Fig. [TJ The sliding s induces the opening 
of the crack, respect to the opening function 9(s). This one depends upon the crack's shape ato (a is a 
physical length). In case of sinusoidal crack shape, naming £ the (metric) abscissa of the crack along 
e 2 , A and p the mathematical amplitude and periodicity, we have (a vanishes in 6): 
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Thermodynamics: state 
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The free energy pty is related to the elastic energy stored by the sound material (with C the Hook 
stiffness tensor). The derivative of it respect to state variables lead to the thermodynamic forces: the 
stress a and the force S (Eq. 0), that has the clear meaning of the projection of the stress o that have 
the will make the crack slide. 

2p^(e,s) = (e-e s ) : C : (e - e s ) (4) 
a = C:e, S = 9'(s)a u + a 12 (5) 

Elasticity domain for a crack 




Figure 2: Yield surface (active crack), for any s and for s = (light gray) 



The yield surface f(a, S) is defined (in the thermodynamic forces space) as a limitation of the 
force S respect to two phenomena: one cohesive and the other frictional. The fist, is a Barenblatt 
adhesive term B ^ that represents cohesion, i.e. when crack does not exist; this cohesion can be 
seen as the cement paste role in a concrete. The second member introduces the term ip ^ that can 
be seen as a Coulomb's friction (at the macro scale level). 

f(a,S) = \S\-B + pa 11 (6) 

The elasticity domain is then defined as / < and the dissipative evolution (siding of the fracture) 
occurs when / = 0. Using Eq. ©, the combination of them, / ^ 0, define the yield function (Eq. [7]) 
and the admissibility condition (Eq. [8]), in the coordinates: 

\9'(s)a u +cr 12 | < B - pa n (7) 



The elasticity domain is represented onto Fig. [21 Non admissible domain is hatched. It can be 
seen as the case when tension becomes large enough in order to separate crack's faces. When B = 0, 
the yield surface tends to a Mohr-Coulomb's one, inclined respect to the slope 9'(s). 

Due to mathematical properties of the crack shape (Eq. [3]), the opening function 9(s) is even and 
the maximum M = max(#'(s)) is obtained for s = kp, then for s = 0. Then the derivative at s = 
belongs to sub-differential defined as 0'(0) G [— M, M\. The corresponding yield surface is non con- 
vex (Fig. [2]); it can be seen as the collection of all the possible yield surface for any s. 



Thermodynamics: dissipation 

For shake of simplicity we retain an associated model: the potential of dissipation is the elasticity 
domain. The flow rule is deduced from the normality rule (in which A ^ 0) z.e. s — f/S X — sign (S) A. 
Then comes that S and s have the same sign, proving that the thermal dissipation Eq. [9] is positive 
under the admissibility condition (Eq. [8]). 

V = Ss=(B -ipa n )\s\ ^0 (9) 

Yield stress in the tension and compression case 



The stress is o = & n <g> n; with n = cos(a)e*i + sin(a)e 2 and (with no restriction) ^ a ^ tt/2. 
In the specimen's orthonormal basis [ft, m], the axial strain writes (from[T][2]): 



a 



£ nn = + Q(s) cos 2 (a) + s sin (a) cos(a) (10) 
E 

Following in some way [CO, we consider that all the cracks are possible in the material; the one 
that will open first will be the one whose orientation a leads to the highest value of the yield function. 
The condition of yielding / = (Eq. [7]) writes: 

|a||0'(s) cos 2 (a) + sin (a) cos(a)| = B — ipacos 2 (a) (11) 

The positiveness of the second member represents the admissibility condition (Eq. [8]); it is always 
true in compression but represents a limit in tension. It can be shown that (Eq. [9]) yields s ^ 0, 
in tension and s ^ in compression. Then, the crack initializes while considering 6 ) '(0 + ) = M in 
tension and #'(Cr) = — M in compression. The crack angle and the stress value is obtained in the two 
cases while minimizing the (absolute) stress level and we obtain (with superscript * for tension and c 
for compression: 



a 1 



2B. k 1 



tan(2a t ) = — — - (12) 



(f + M+ y/l + {(p + M) 2 ip + M 

^ = 2 B tan(2« c ) = (13) 

ip + M- y/l + (if + M) 2 V + M 

The ratio — a c /a 1 appears to depend only upon the value of ip + M. Considering a Coulomb's 
friction cp = 0.5 and M — 1 (approximatively relevant of measured roughness jH), this ratio is 10.9, 
typical for concrete H. The angle of cracks are respectively a 1 =17 and a c =73 degrees for tension 
and compression (Fig. 4) and these value are reasonable respect to experiments. 



Evolution 



The post-peak response can then be obtained from (Eq. [TTj) . as angles a are set. The Barenblatt 
constant B is supposed to depend upon the sliding s, in order to represent the damage of the material. 
Using a linear dependancy, we suppose that no more cohesion remains after an half step: 

o A 2max t | g |\ + 



The obtained stress to strain curve (for p = 5 10~ 3 , E = 70 GPa, B = 6 GPa, Fig. 3) exhibits 
softening in both tension and compression. The zero stress state is obtained for a larger strain in 
tension than in compression. These two facts lead to consider that the proposed model is only relevant 
to yield stress but not yet to describe anelastic behavior. 




Figure 3: Tension to compression curve Figure 4: Crack kinematics 

Conclusion 

This approach seems to be able to predict yield stress and crack angle in tension and compression. 
An extension in case of bi-tension or compression will be done. The anelastic behavior needs to 
be enhanced in order to be predictive. Softening behavior will have to be considered with respect 
to some material scale in order to avoid localization problems (mesh dependancy); this micro-scale 
being already present in the model allow to think that it is possible. 
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